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Abstract 

For the multisensor autoregressive moving average (ARMA) 
signals with unknown model parameters and noise 
variances, using the recursive instrumental variable (RIV) 
algorithm, the correlation method and the Gevers-Wouters 
algorithm with a dead band, the fused estimators of 
unknown model parameters and noise variances can be 
obtained. Then by substituting them into optimal fusion 
signal smoother weighted by scalars, a self-tuning 
distributed fusion Kalman smoother is presented. Using the 
dynamic error system analysis (DESA) method, it is 
rigorously proved that the self-tuning fused Kalman signal 
smoother converges to the optimal fused Kalman signal 
smoother, so that it has asymptotic optimality. A simulation 
example shows its effectiveness. 
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Introduction 

Information fusion techniques combine data from 
multi-sensors and related information to obtain more 
specific inferences than those obtained by using a 
single sensor [Liggins et al. (2009)]. The existing 
information fusion Kalman filtering is mainly focused 
on the information fusion Kalman filtering with 
known model parameter and noise statistics. However, 
in many applications, the model parameters and/or 
noise variances are usually unknown. The filtering for 
the systems with unknown model parameters and/or 
noise variances is called self-tuning filtering [Deng 
( 2012 )]. 

Recently, several self-tuning weighted fusion filters 
[Deng et al. (2008), Deng & Li (2007), Sun (2007)] have 
been presented, whose drawbacks are that only the 
noise variances are assumed to be unknown, while the 
model parameters are assumed to be known. In the 
existing results, most of them are for the single- 
channel signal systems with white measurement noise 
[Tao et al. (2010), Liu & Deng(2010), Gao & Deng 


(2010)]. And, up to the present, the information fusion 
is mainly focused on filtering fusion for the single- 
channel multisensor systems with white measurement 
noises, while smoothers for ARMA signals are rarely 
reported[Tao & Deng (2010), Tao & Deng (2012), Deng 
& Li (2007)]. 

In this paper, for multisensor autoregressive moving 
average (ARMA) signals with a common disturbance 
noise, using the Kalman filtering method, a self-tuning 
distributed information fusion Kalman smoother 
weighted by scalars is presented for the case when 
model parameters and noise variances are both 
unknown. By the dynamic error system analysis 
(DESA) method [Ran et al. (2009)], it is proved that the 
self-tuning distributed fusion Kalman smoother 
converges to the optimal distributed fusion Kalman 
smoother with probability one, so it has asymptotic 
optimality. Compared with[Tao & Deng (2010), Tao & 
Deng (2012), Deng & Li (2007)], the paper uses 
different signals model with a common disturbance 
noise. In addition, unknown model paremeters and 
noise variances are different from [Tao & Deng (2010), 
Tao & Deng (2012), Deng & Li (2007)]. 

Problem Formulation 

Consider the multisensor single-channel ARMA signal 
with L -sensor 

A(cf')s(t) = C(cf')w(t) ^ 

y i (t) = s(t) + v i (t)' i = ( 2 ) 

vft) = f(t) + eft) 

where t is the discrete time, s(t) is the ARMA signal to 
be estimated, yft) is the measurements of the zth sensor, 
w( t ) and vfi) are the input noise and measurement 
noises of the/th subsystem, respectively, v i (t) contain a 
common disturbance noise £(t) , q~ 1 is the backward 
shift operator, cf l s{t) = s{t- 1) , A{q~') and C(q~') are 
stable polynomials having the form as A(q ') = 1 + 
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a x q ' +---a„ a q and C(q ') = \ + c x q 1 t — c n q " c ,n a is 
the order of polyno-mial A (q ' ) , n c is the order of 
polynomial of C{q~'). 


Assumption 1 w(t) , g(t) and e,(t) (i = 1,2, •••,£) are 
uncorrelated white noises with zero mean and 
variances and a 2 , respectively, and the input 

noise w(t) and the measurement noise v i (t) have the 
relation 


E- 




o - ,; 

o 


aj+trlS,, 


\S, t 


(4) 


where the symbol E denotes the mathematical 
expectation, the superscript T denotes the transpose, 
S tt =l,S lk =0(t*k). 


From (3), it can be seen that the measurement noises 
v i (t)(i = 1 are correlated,i.e. 

R,, = E[v,(0vJ(0] = o-f(i^./),R, =E[v,(0v/(0] = cr^ + cr (5) 

Assumption 2 A(q ')and C(q 1 ) are stable polynomials. 


Assumption 3 The model parameters a i , c , , and noise 
variances cr 2 , cr; and cr 2 are unknown. 

Assumption 4 The realizations of measurement 
stochastic process y ; (7) (z = 1,2, — ,Z,) are bounded for t , 
with probability one. 


The Optimal Fusion Kalman Signal Smoother 
Weighted by Scalars 

Setting w(t) = iv(t-l) , yields that w(t) has the variance 

The ARM A signal (1) can be rewritten as 

A(q x )s(t) = C(q-')w(t) (6) 

C (q ) = c x q + • • • + c n _q ~"< , c, = c w , / = 1,2, , 

"c=n c + 1 ( 7 ) 

From (2) and (5), we have the state space model with 


the companion form 

x(t + 1) = 0x(t) + Iw(t) (8) 

y i {t) = Hx{t) + v i (t), i = lX-,L (9) 

s(t) = Hx{t) (10) 

-a x 

0= '■ h a -\ // = [l 0 ■■■ 0] 

~a„ a 0 O' 

^M 1 C 1 ”• c n a \, C J = 0(i > n c ) (11) 


Lemma 1 [Deng (2012)]. For the multisensor system 
(7)-(9) with known model parameters and noise 
variances, the zth sensor subsystem has the local 
steady-state optimal Kalman predictor 

*i(t + 1 1 1) = I f / Pi (t)x i (t \t-\) + K Pi {t)y i (t) (12) 

P Pl (t) = 0-K Pi (t)H (13 ) 

K Pi (t) = 0S l (t\t-\)H J Qj(t) (14) 

Q a (t) = HE l (t\t-\)H 1 + R t (15) 

where the local prediction error variance matrix 
E t (t 1 1 - 1) satisfies the Riccati equation 


Z i {t + \\t) = 0[Z i (t\t-\)-Z i (t\t-\)H 1 (HZ i (t\t-\)x 

H t +R i y > HI i (t 1 1 - l)]0 T +ra;,r T (i 6) 

and the local prediction cross-covariance 

I lj (t\t -\) = E[x i (t\t -X)^ (t\t -\)}, 
with x t (t | r — 1) = x(t) — x x (t 1 1 — 1) satisfies the Lyapunov 
equation 

z,j(t+\\t) = r Pi (t)z y (t\t- w T Pi (o + ra 2 w r T + 

K pi (t)RjK T Pj (t ) i±j i,j = \,2,-,L (17) 

/ / 

Lemma 2 [Deng (2012)]. For the multisensor system (1) 
and (2) with known model parameters and noise 
variances, the ^ sensor subsystem has the local 
optimal Kalman signal smoother as 

N 

x i (t-N\t) = x i (t~N\t-N-\) + 'Y J K i (t-N\t-N + j)x 

i = 0 

£ i(t ~ N + j),i = (18) 

The optimal Kalman signal smoother is given by 

s i (t-N\t) = Hx i (t-N\t) (19) 

where the time-varying optimal Kalman predictor 

x t (t - N 1 1 - N-l) can be computed by (12), and we 
have 

s i {t) = y i {t)-Hx i {t\t-\) (20) 

i - 1 

K,(t\t + j) = Z, (t 1 1 - - 1) {]] Wl (t + k)}H T Q-J (t + j) (21) 

k = 0 

and the cross-covariance among smoother errors 
s i (t-N\t) = s(t-N)-s i (t-N\t) / is P s ..(t-N\t) = HP ij {t-N\t)H T , 
and the error variance matrices and covariance 
matrices of local Kalman smoother are given as 

N 

P i (t-N\t) = E i (t-N\t-N-l)-^K i (t-N\t-N+j)x 

7=0 

Q a (t-N + j)Kj(t-N\t-N + j) (22) 


P..(t -N\t) = P iN {t -N)E y (t-N\t-N - 1 )P] N {t -N) + 


v ' 

N 


yk(h^) KUt-N) 


L ip 

P = 0 


0 
( b 

KJ p (t - N) 

1 

0 

1 

^ 1 


(23) 
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with the definition P u (t - N \ t) = P t (t — N \ t) . 

N 

'F m (t-N) = I n -Y,K i (t-N\t-N + k)Hx 

k = 0 

P p ,(t-N + k,t-N). 

/ 

N 

K;(t-N) = ~Y J K i (t-N\t-N + k)Hx 

k-p+\ 

W pj (t-N + k,t-N + p + l)r p = 0,---,N-l; 

/ 

N 

KJ p {t-N) = YjK^-Nlt-N + kWx 

k=p + 1 

'F pi {t-N + k,t-N + p + \)K pi (t-N + p)- 

K t (t-N\t-N + p) p = 0,—,N-\\ 

/ 

K”(t-N) = O r K v „(t-N) = -K i (t-N\t) (24) 
where ¥G (t + k, t) = 'F pi (t + k- 1 ) • • • P pi (. t ), <F pi (t,t) = /„ . 

Lemma3 [Deng (2012)]. For the multisensor system (7)- 
(9) with the assumptions 1 and 2, the optimal 
information fusion Kalman signal smoother 
s 0 (t-N\t) is given by 

L 

s 0 (t -N\t) = '^(O i (t-N\ t)$i(t -N\t) 

where the optimal scalars weighting coefficient vectors 
C 0 i(t-N\t) i s given by 

Wt-N\t), -, co L {t-N\t)}=[e\P;\t-N\t)Y l x 

e]- l e T (P s (t-N\t)y > (26) 

where e T = [1, — ,1] , defining the LxL matrix as 

P s {t-N\t) = (tr P sij (t-N | 0 ) , ij = 1,2,-, L (27) 

Self-tuning Fusion Kalman Smoother Weighted 
by Scalars 

When model parameters and noise variances are 
unknown, substituting estimators into the optimal 
fusion Kalman smoother will yield the self-tuning 
fusion Kalman smoother which consists of the 
following steps: 

Step 1. Applying the system identification algorithm 

[Ljung (1999)], the estimator a,(l) can be obtained. 

Applying the correlated method [Ran et al. (2009)] and 

Gevers-Wouters algorithm [Gevers & Wouters (1978)], 

fusion estimators c,.(7), ay.L) dl{t) d 2 (t) R ,,(0 and 

/ / / 

R, (t) can be obtained, where R tj (t) = dl{t) {i ^ j) , 
R i (t) = dl(t) + d 2 (t) . And the model parameters and 
noise variance estimators are consistent, i.e. 

a,(t) -> a,., c,.(f) -> c,. al„d 2 Jt) -> o'! 


(!) -> a\ R, (0 -> /?. R,y (0 -> ^ 

/ / 

as f -» °°, w.p.l (28) 

Step 2. Substituting the estimators d, (i) and c,(7) into 

(11) yields <Z>(r) and /’(r) . In Lemma 1,0 r o 2 v R, 

/ / / 

and R, y (i,j = !,■■■, L) are replaced by 0(t) k(t) <?l(t) 

/ / / 

R,(?) and R (t) , respectively. And the estimators 
Z i (t 1 1 -1) satisfy the self-tuning Riccati equations 

i; (t + 1 1 o = Spaz (t\t-\)-z(t\t - i )H J (H x 
4 (t 1 1 - 1)// T + R, (*))' 1 HZ, (t\t- 1)]0 J (t) + 

(29) 

and the prediction cross-covariance matrices satisfy 
the self-tuning Lyapunov equations 

i,(i+i 1 1) = 5^(03, (* 1 1 - 1 )^( 0 + ht)6-l(t)f\t)+ 


K^RjWlit), i*j,i,j = ia,-,L (30) 

with the definition Z u (t + 1 1 1) = Z t (t + 1 1 1) . Substituting 
the estimators into (12), the self-tuning local Kalman 
predictor can be given as 

x- (t + \\t) = P R , (t)xl (t\t-\) + K Pi ( t)y it) (31) 

P R ,(t) = 0(t)-k R ,(t)H (32) 

k Pi (t) = 0(t)Z i (t\t-X)H 1 Qj(t) ( 33 ) 

Q a (t) = HZ(t\t-])H J + R i (t) ( 34 ) 


Substituting the estimators into (18)-(24), yields the 
self-tuning local Kalman smoother x-(t-N\t), 
( i = 1 ,---,L) as 

N 

x s i (t-N\t) = x*(t-N\t-N-l) + Yk i (t-N\t-N+ j)x 

j = 0 

£ i (t~N + j),i = !,■■■, L ( 35 ) 

£i(t) = y i (t)-Hx s i (t\t- 1 ) ( 36 ) 

i - 1 

k, ( t\t + j) = Z, hlt-mn ^ ( t + k)}H T (Hx 

k = 0 

Z,(t + j\t + j-\)H'+R,(t))-' ( 37 ) 


The error variance matrices and covariance matrices of 
self-tuning local optimal Kalman smoother are given 
as 




j = o 


Q,(t-N + j)Kj(t-N\t-N + j) 


(38) 


PM-N\ t ) = <P m (t-N)ZM-N\t-N-W] N {t-N) + 
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~KU) o ' 

1 

1 

0 R iJ (t)_ 

kf p (t-N) \ 


J^KZit-N) K; p (t-N) 

P = 0 

with the definition P..(t -N\t) = P i (t-N\t). 

N 

'r m (t-N) = I n -'££ i (t-N\t-N + k)Hx 


(39) 


' F pi (t-N + k,t-N ). 

/ 

N 

K”(t-N) = -^K i (t~N\t-N + k)Hx 

k-p+l 

'F pi (t-N + k,t-N + p + \)f(t) p = 0,-,N-V, 

/ 

N 

K] p (t — N) = ^ K i (t-N\t-N + k)Hx 

k=p + 1 

^.(t-zv+^t-zv+p+D^.^-zv + p)- 

K i (t-N\t-N+ p) p = 0,---,ZV-l; 

/ 

k>-zv) = o <;(t-zv) = -l,.(t-zv|0 
/ 

where (/ + *,#) = ^ (0, ^ (Z, 0 - 4 - 


And the self-tuning local Kalman signal smoother is 
given by 

s s l (t-N\t) = Hil(t-N\t) ( 41 ) 

Step 3. Applying (26) and (27), the estimates 
P s (t-N\t) and fit, (Z - ZV 1 1) are obtained, and the self- 
tuning fused Kalman smoother is given as 

L 

K (t - N 1 1) = ^ fit,. (t-N\ z)y s (t-N\t) ( 42 ) 

1=1 V 

The above three steps are repeated at each time t. 


Convergence Analysis 

Convergence of Self-tuning Lyapunov equation 
Theorem 1. For the multisensor system (1) and (2) with 
the assumptions 1-4, the solution -^(Z + l I Z) of the 
self-tuning Lyapunov equation (30) converges to the 
solution 27,./ (t + 1 1 1) of the optimal Riccati equation (17), 
i.e. 

+ 1 1 t)-Sy(t + \ | Z)] -» 0, as ? oo,i.a.r (43) 

Proof. It has been proved that the solution 27,.(Z + 1|Z) 
of the self-tuning Riccati equation converges to the 
solution 27,. of the steady-state optimal Riccati 
equation [Tao & Deng (2010)], i.e. 

2/ (Z + 1 1 Z) - 27 — > 0 , as t — > oo , i.a.r (44) 

According to the stability theory of the Kalman 
filtering [Kamen & Su (1999)], the solution 27,.(Z + 1|Z) 
of the optimal Riccati equation converges to the 


solution 27,. of the steady-state optimal Riccati 
equation , i.e. 

I ft + 1 1 Z) - 27 -» 0 ( 45 ) 

From (44) and (45), we have that 

4(Z + 1 1 Z) - ^.(Z + 1 1 Z) = (27, (Z + 1 1 Z) - 27,.)- 
(27,.(Z + 1 1 Z)-27,.) — > 0, as Z — >■ oo, i.a.r (46) 

From (14) and (33), according to the stability theory of 
the Kalman filtering [Tao & Deng (2010)], the steady- 
state optimal Kalman predictor gain matrix K Pj i s a 
continuous function with respect to cp , 27. and , i.e. 
K Pi = OZ i H r {HZ i H 1 + R i )~ l =f i {0,E i ,R i ) . From (11) 
and (28), we have <P(Z) — > &,f(t) — > r. According to 
the continuity of f , applying (28) and (46) yields that 

{K Pi (t) - K pi {t)] — > 0, as Z ^oo, i.a.r (47) 

Applying (13), the steady-state predictor transition 
matrix f'p, is a continuous function with respect to K pj 
and (p , i.e. f'p, = & ~ K Pi H = gfK Pi ,&) . From (47) and 
the continuity of g, , we can yield that 

[E Pl (t)-P Pl (t)]^0, as t ->oo, i.a.r (48) 

Setting P n (t) = + A E P ,(t ) , k„(t) = K„(t ) + dK„(Z) , 

then 

A 'Fpft) — » 0, dK p ,-(t) — > 0 , as t —> oo , i.a.r (49) 
Subtracting (17) from (30), and setting 

£,(Z)=i,(Z + l|z)-27,(z + l|0 ( 50 ) 

we obtain the dynamic variance error Lyapunov 
equation 

Eyf) = 'F Pi {t)E ij (t-\)'Ff j (t) + U ij (t) (51) 

with the input 

u t] {t) = A Ep.m.ft 1 1 - Wfft) + A 'FnWjit | Z - 1) x 
A Ppjit) + n,(z)T iy 0 1 1 -1)A $1 it) + - (52) 

Col,r J +k P ft)Ryt)r Pj {t)-K P ft)R ij Klit) 

According to the stability theory of the Kalman filte- 
ring[Kamen & Su (1999)], applying Lemma 6 yields 
that K Pj {t), 'Epft), l P Pi {t) and K p ft) are uniformly 
asymptotically stable. Applying the boundedness of 
r(t) , 6- 2 w {t) yields that r(t)olft)r T (t) is bounded. 
Similarly, it can be obtained that Kp,.(t)R, y (t)Kp ; (t) is 
bounded. Applying Lemma 3 [Tao & Deng (2012)] to 
(30) yields 27. (z | Z — 1) is bounded. Therefore, from (52), 
applying (28), (43) and (49) f(t) — > r yields that 
Uy{t) — > 0, as Z -> oo , i.a.r . Applying Lemma3 [Tao & 
Deng (2012)] to (51) yields 
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Ejj(t) — > 0, as t— »oo, i.a.r (53) 

which gives (43). The proof is completed. 

Convergence of the Self-timing Local And Fused Kalman 
Smoothers 

Theorem 2. For the multisensor systems (1) and (2) 
with the assumptions 1-4, the local self-tuning Kalman 
predictor x-(t + 1 1 1) converges to the local time-varying 
optimal Kalman predictor x i ( t + 1 1 1) in the sense that 

[y s (l + l|l)-x,.(i + l|l)]— ) >0, as ? “ ^ °°, i.a.r. (54) 

Proof. Denoting 5.h) = x/(t + 1 \t)-xft + \ \ t), subtracting 
(12) from (31) yields a dynamic error system 

8ft) = V'pi m (t - 1) + w, (0 (55) 

u, (t) = A W Pi (t)x, (t 1 1 - 1) + AK Pi (t)y, (. t ) (56) 

From the boundedness of K p ft ) and the assumptions 
4, K p ft)y ft) is bounded. Since P P ft) is uniformly 
asymptotically stable, hence, the application of Lemma 
4 [Tao & Deng (2012)] to (31) yields that x j {t\t-L) is 
bounded. From the boundedness of x. (t\t-\) and 
y, (0 , applying (49) to (56) yields that 

uft) -» 0 , as t — > oo , i.a.r. (57) 

Since Ppft) is uniformly asymptotically stable, 
applying Lemma 4 [Tao & Deng (2012)] to (55) yields 
that 

8ft) -» 0 , as t -> oo , i.a.r. (58) 

Hence, (54) holds. The proof is completed. 

Theorem 3. For the multisensor system (1) and (2) with 
the assumptions 1-4, the self-tuning fused Kalman 
smoother y, (/ - N \ I ) converges to the optimal fused 
Kalman smoother x 0 (t ~ N \ t) | n a realization, i.e. 

[jcJ(t-tV|f)-jc 0 (f-tV|t)]— ; >0, as t— >oo, i.a.r. (59) 
Proof. From (20) and (36), applying (54) yields 

[sft)-£ft)\— >0, as t— ><», i.a.r (60) 

From (21) and (37), applying (28), (43) and (48) and the 
property of the continuous function yields 

[kft-N\t-N + j)-Kft-N\t-N + j)-\^ 0, 

as t — » oo , i.a.r (61) 

Setting s, (t) = e, ft) + As, (t ) , Kft - N \t - N + j) = 

Kft — N\t — N + j) + AKft — N\t — N + j), from (60) and 
(61), we have 

As ft) -» 0, AKft -N\t-N + j)->0 

/ 

as t — > oo , i.a.r (62) 

Denohng d t (t) = x- (t - N \ t) - x, ( t-N\t ) , and subtractting 


(18) from (35) yields 

d, (t) = x- (t - N \t - N - 1) - x, (t - N \t - N - 1) + ?/, (t) (63) 

N 

1ft) = ^ft — N\t — N + j)sft -N + j)- 

7=0 

Kft-N\t-N + j)sft-N + j) (64) 

From (64) we have 

N 

1 ft) = 'Y J Kf t ~ N \t~ N + j)Asft -N + j) + 

7=0 

AKft-N\t-N + j)sft-N+j) + 
AKft-N\t-N + j)Asft-N + j) (65) 

From (36), applying the boundedness of yft) and 
x) (/ 1 / - 1 ) yields that eft) is botmded. From (65), 
applying (62) yields 

lift) -» 0 , as t -> oo , i.a.r (66) 

From (54) and (65), we have 

dft) — > 0, as t — » oo , i.a.r (67) 

From (22), (23), (38),(39), applying (28), (46)-(48) and 
(61), we have 

Pft-N\t)^Pft-N\t) P j (t-N\t)^P ij (t-N\t) 

/ 

as t — >oo, i.a.r (68) 

From (26) and (27), each element of a>ft-N\t) j s a 
continuous function of elements of P kJ (t-N\t) r we 
have 

tbft- N\t)-> (oft -N\t), as t -> oo, i.a.r (69) 

Setting coft-N\t) = coft-N\t) + Aa> ft- N\t) , from (69) 
we have that Ad) ft - N\ t) — > 0 , subtracting (25) from 
(42) yields 

L 

x'o(t — N 1 1) — x 0 (t — N |t) = ^ ®,. (t - N 1 1) x 

j=i 

L 

[x-(t -N\t)-xft-N\t)] + Y^ A ®ft - N 1 C t-N\t ) (70) 

i=l 

Applying (69), Ad> ft - N \ t) — > 0 , and boundedness of 
Xj (t-N\t) yields that (59) holds. 

Simulation Example 

Consider the ARMA signal s(t) with measure-ment 
noises 

A{q- ] )s(t) = C(q-')w(t) (71) 

yft) = s(t) + f(t) + eft), i = 1,2, ■■■,L (72) 

where f(t) , eft) and w(t) are independent Gaussian 
white noises with zero mean and variances , a\,. 
and cr^, . In simulation, we take 

A(q 1 ) = l + ^j^ 1 + a^q ” = 1 + 0.2 q — 0.48^/ ”, 
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C(q- l ) = l + Cl q- 1 =l + 0.4q-\al =0.1 r al = 0.7, 

^ = 0.1, <4 = 0.2, <4 = 0.3. 

The problem is to design the self-tuning information 
fusion filter of s(t) ,where ay , cr: t a^(i = 1,2,3) , a, , a 2 
and c, are unknown. 

The simulation results are shown in Figure 1-Figure 6. 
From Figure 1-Figure 5 in which the curved lines 
denote the estimates, and the straight lines denote the 
true values, it can be observed that the estimates of the 
model parameters and noise variances are consistent. 
Figure 6 shows that the self-tuning fusion Kalman 
signal smoother s*(f-l|f) converges to the optimal 
fusion Kalman signal smoother s 0 (f-l| t ) , i.e. its error 
converges to zero. 
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FIG. 5 THE FUSED ESTIMATE CURVES OF MODEL PARAMETER 
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FIG. 3 THE FUSED ESTIMATE CURVE OF NOISE VARIANCE 


FIG. 6 THE ERROR CURVE BETWEEN THE SELF-TUNING AND 
OPTIMAL FUSION KALMAN SIGNAL SMOOTHERS 
e(t) = s* 0 (t-l\t)-s 0 (t-l\t) 

Conclusion 

For multisensor ARMA signal with unknown model 
parameters and noise variances, by the system 
identification, correlation method and the Gevers- 
Wouters algorithm with a dead band, a self-tuning 
distributed information fusion ARMA signal Kalman 
smoother was presented. By the DESA method, it was 
rigorously proved that it converges to the optimal 
fusion ARMA signal Kalman smoother with 
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probability one, so that it has asymptotic optimality. 
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